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摘要 ” 喷 淋 雾 化 被 广泛 应 用 于 核 动 力 设备 中 ， 以 保证 其 安全 稳定 运行 。 喷 淋 时 ， 部 分 液 滴 会 撞击 设备 壁面 ， 
以 液 膜 的 形式 进行 扩展 和 演变 。 基 于 Fluent DPM-to-VOF 进行 了 液 滴 射 流 撞 壁 研究 ， 使 用 两 液 滴 滴 落 过 程 实 
验 图 像 验 证 了 模拟 方法 的 准确 性 ， 研 究 了 液 滴 射流 注入 速度 、 入 射 角度 对 壁面 液 膜 形态 的 影响 规律 ， 得 到 了 
液 滴 溅 射 率 随 入 射 条 件 的 变化 情况 ， 并 分 析 了 液 滴 撞 击 壁 面 时 重力 、 表 面 张力 、 速 度 分 量 等 对 液 膜 形态 扩展 
过 程 、 液 滴 溅 射 率 的 影响 。 研 究 发 现 ， 增 大 注入 速度 能 加 大 对 液 膜 的 扰动 ， 液 膜 与 壁面 的 接触 面积 显著 增 大 ， 
液 膜 区 宽度 和 高 度 增 加 ; 增 大 入 射 角度 对 液 膜 形态 的 改变 影响 较 大 ， 体 现在 液 膜 区 宽度 增加 和 液 膜 区 高 度 减 
小 ， 液 膜 表 面 稳定 性 减低 。 溅 射 率 随 着 注入 速度 和 入 射 角度 的 增 大 而 增 大 ， 其 中 ， 入 射 角度 对 溅 射 率 的 影 上 
更 大 。 
关键 词 DPM-TO-VOF; 液 滴 射 流 撞 壁 ， 液 膜 扩 展演 变 ， 液 滴 滴 落 融合 
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Study on the Droplet Jet Impingement Process on a Wall Based on DPM-to-VOF 


CHEN Qingshan1， WANG Mingjun!, TIAN Wenxil ， QIU Suizheng!, SU Guanghuil 


1(Nuclear Thermal-hydraulic Laboratory, School of Nuclear Science and Technology, Xi'an Jiaotong University, Xi 'an Shaanxi, China) 


Abstract: [Background]: Spray atomization is widely used in nuclear power equipment to ensure its safe and stable 
operation. [Purpose]: This study aims to obtain the characteristics of the extended evolution of a liquid film formed 
by a droplet jet impacting a vertical wall. [Methods]: This paper conducted a study on droplet jet impinging on a wall 
based on Fluent DPM-to-VOF, and validated the accuracy of the simulation method using experimental images of 
two droplet falling processes. The research investigated the influence of droplet injection velocity and incident angle 
on the wall liquid film morphology, and obtained the variation of droplet splash rate with incident conditions. 
Furthermore, the effects of gravity, surface tension, velocity components, and other factors on the expansion process 
of liquid film morphology and droplet splash rate during droplet impact on the wall were analyzed. [Results]: The 
study found that increasing the injection velocity could enhance the disturbance to the liquid film, significantly 
increasing the contact area between the liquid film and the wall, as well as increasing the width and height of the 
liquid film area. Increasing the incident angle had a significant impact on the change of liquid film morphology, 
reflected in the increase in the width of the liquid film area and the decrease in the height of the liquid film area, 
leading to decreased stability of the liquid film surface. The splash rate increased with the increase of injection 
velocity and incident angle, with the incident angle having a greater impact on the splash rate. [Conclusions]: The 
study of jet impingement on the wall can lay the foundation for the subsequent study of wall cooling under the 
spraying. 
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雾 化 和 喷 淋 在 了 


作用 。 


和 壁 


喷 淋 系统 的 使 


制 放 射 性 排放 等 多 个 功能 。 
膜 的 情况 ， 


系统 


液 滴 撞 


Model, DPM) 到 连续 相 (Volume of Fluid, VOP) 的 耦合 模型 
面 、 壁 膜 的 过 程 。 

为 了 实现 离散 相 到 连续 相 的 耦合 ， 国 内 外 学 者 进行 了 
实验 平台 进行 了 多 液 滴 滴 落 聚 并 实验 下 
并 预测 了 线 和 膜 保 持 连续 、 
-DPM-VOF 集成 算法 来 求解 结 冰 过 程 的 相间 传 热 问题 ， 
结 冰 过 程 
分 析 ] 


等 人 利用 液 滴 聚 # 
续 相 液 线 和 液 膜 的 过 程 ， 


ja 


ri 


提出 了 一 


性 和 可 靠 性 ， 揭 示 
拟 了 静止 液 滴 的 融合 过 程 ， 
拟 运 动 液 滴 的 碰撞 蜗 


[ 业 、 农 业 、 消 防 和 
j 对 核 动力 设备 的 
忌 喷 淋 、 稳 压 器 喷 淋 和 安全 壳 喷 淋 在 内 的 多 个 喷 淋 系统 ， 起 到 了 排出 堆 世 余热、 冷凝 燕 汽 、 降 温 
在 这 些 喷 淋 的 实际 应 用 场景 中 ， 有 时 会 出 现 液 相 以 液 滴 的 形式 撞击 壁面 、 
典型 的 例子 如 稳 压 器 喷 淋 系统 中 液 滴 撞 击 到 稳 压 器 内 壁面 和 稳 压 器 液 相 空间 、 安 全 壳 喷 淋 
击 堆 内 设备 的 壁面 等 。 对 于 此 类 应 


医疗 等 领域 应 用 广泛 ， 特 别 是 在 核电 广 和 核 动力 设备 中 发 挥 着 重要 
安全 运行 和 反应 堆 的 完整 性 有 着 十 分 重要 的 作用 凹 ， 设 置 了 包括 堆 
降 压 和 限 
液 面 


的 模拟 研究 ， 可 以 使 用 一 种 由 离散 相 (Discrete Particle 


昌 ， 从 而 精准 捕捉 到 离散 液 滴 进入 液 面 或 者 撞击 壁 


广泛 研究 ， 并 进行 了 大 量 的 试验 和 模拟 。DaliliDl 
究 ， 观 察 到 粘性 液 滴 在 固定 钢板 上 的 聚 并 形成 连 
发 生 断 裂 的 条 件 。LiuBl] 等 人 提出 了 一 种 固定 网 格 烩 
并 通过 两 个 甘油 液 滴 滴 落 案例 验证 了 耦合 方法 的 准 
液 滴 的 形成 、 运 动 以 及 相 变 、 传 热 的 细节 。 刘 栋 外 等 人 采用 FTM 方法 模 


粘性 、 液 滴 尺 寸 和 环境 因素 对 融合 过 程 的 影响 ， 并 发 展 了 SPH 模型 来 模 


合 过 程 ， 


再 现 了 独立 液 滴 到 连续 液 滴 的 融合 过 程 。 Zhang51 等 人 利用 CONVERGE 代码 


' 拉 格 朗 日 - 欧 拉 混合 过 渡 模 弄 


de 


苗 述 了 喷 淋 液 滴 冷 却 成 膜 的 全 部 过 程 , 利用 从 多 个 液 滴 滴 落 过 程 


型 (TFM) 的 详 


生 连 贯 水 膜 时 


的 实验 数据 对 模型 进行 了 验证 ， 此 外 ， 还 对 喷雾 撞击 形成 液 膜 的 过 程 进 行 了 过 渡 模 型 


比较 。 而 


和 基于 颗粒 的 薄膜 模 


生 的 铺展 形态 ， 唐 亮 中 等 人 利用 实验 研究 了 倾斜 射流 撞击 壁面 形成 液 膜 的 关键 特 
等 人 研究 了 不 同 壁面 对 汽油 液 滴 撞 壁 结果 转换 的 影响 。 


完 人 员 同 样 进 行 了 大 量 的 射流 撞 


完了 射流 撞击 玻璃 表面 产 
征 ; 及 振 山 的 


i 


究 。WangIg 等 人 下 


在 稳 压 器 喷 淋 过 程 中 ， 部 分 液 滴 冲击 稳 压 器 壁面 ， 并 逐渐 聚集 、 铺 展 成 连续 的 液 膜 ， 随 后 在 下 游 区 域 


收缩 、 成 股 的 下 流 。 
拟 方法 对 液 滴 射 流 撞 避 
模拟 计算 ， 并 通过 两 液 滴 滴 落实 验 的 详细 情况 对 模型 进行 了 模型 验证 。 


研究 方法 和 理论 模型 
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1.1 DTV 方法 


对 于 射流 液 滴 撞击 壁面 和 
的 耦合 。DTYV 方法 在 以 下 两 种 情形 发 生 转 换 : 


首先 ， 当 离散 颗粒 撞 避 


液 膜 的 形态 和 铺展 过 程 都 会 影响 其 与 稳 压 器 壁面 的 传 热 特性 ， 因 此 采用 合适 的 数值 模 
Ff 壁面 的 模拟 很 有 必要 。 本 文 基于 DPM-to-VOF (DTV) 方 法 进行 了 多 液 滴 撞击 壁面 的 


的 研究 方法 ，Fluent 提供 了 DTYV 方法 ， 可 以 实现 由 DPM 到 VOF 的 过 渡 过 程 


转换 为 VOF 液 相 。 


其 次 ， 当 离散 颗粒 接触 到 


中 进行 转换 。 


二 


~ 


干燥 壁面 时 ， 会 形成 拉 格 朗 日 壁 膜 ， 当 膜 高 度 超过 设 定 的 闵 值 时 ，DPM 壁 膜 会 


用 户 可 以 指定 DPM 到 VOF 的 过 渡 准 则 ， 其 
等 。 如 图 1 所 示 ，DPM 粒子 与 VOF 液 芯 接触 (a)， 在 计算 网 格 内 ， 接 


欧 拉 壁 膜 的 气 液 界 面 时 ， 若 满足 设 定 的 转换 标准 ，DPM 会 融入 到 VOF 液 相 


包括 VOF 的 相 体积 分 数 阔 值 和 拉 格 朗 日 壁 膜 转 换 参 数 
由 的 VOF 相 体 积分 数 超过 设 定 的 阔 


时 ， 粒 子 的 轨迹 计算 将 停止 ， 
望 将 拉 格 朗 日 


壁 膜 也 转换 为 VOF， 则 可 以 设 定 一 个 膜 高 度 上 的 闵 值 。 


1.2 ”理论 模型 


并 转换 为 VOF (c)， 并 同时 进行 粒子 的 质量 和 动量 的 传递 。 此 外 ， 如 果 和 希 


壁面 上 的 


欧 拉 壁 膜 的 质量 守恒 方程 为 ; 
9 p h E 
人 .(p 1 hV 1 )=m 车 #(1) 


式 中 ，P ,为 液体 密度 ，h 为 液 膜 高 度 ，V ”， 为 平均 液 膜 速度 ，m”。 为 壁面 


肖 
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壁 膜 的 动量 守恒 方程 为 
: +Vs (ph +D vy )=-hnv, P 1, +p 1 hg ， +5 


2 


其 中 ， 式 (2) 左 侧 分 别 为 液 膜 的 瞬 态 效应 和 对 流 效 应 ， 右 侧 第 一 项 表示 气体 压力 、 重 力 垂直 于 壁面 的 分 量 、 
表面 张力 的 影响 效应 ， 第 二 项 为 重力 在 平行 于 壁面 的 分 量 的 影响 效应 ， 第 三 项 和 第 四 项 为 粘性 剪 切 力 影响 
效应 ， 第 五 项 为 液 滴 源 项 ， 第 六 项 为 表面 力 效应 。 

欧 拉 壁 膜 模 型 可 以 通过 改变 膜 方程 的 源 项 来 实现 DPM 和 VOF 多 相 流 模型 的 相互 作用 。 当 离散 粒子 碰 
撞 壁 膜 被 吸收 时 ， 粒 子 的 质量 和 动量 会 分 别 添加 式 (1)、(2) 的 源 项 中 。 


质量 源 项 为 : 
m 8 =m p #(3) 
其 中 ，m ”,。 为 离散 粒子 撞击 液 膜 的 流量 。 
动量 源 项 为 : 
q 和 =m , (V ， Ss ) #(4) 


其 中 ，V。 。 为 离散 粒子 的 速度 。 

当 DPM 离散 粒子 在 撞击 壁面 时 ， 若 壁面 上 不 存在 VOF 连续 相 ， 粒 子 也 将 被 收集 形成 拉 格 朗 日 壁 膜 。 
当 拉 格 朗 日 壁 膜 聚集 并 且 厚 度 超过 设 定 的 最 大 阔 值 时 ，DTV 模型 会 使 之 发 生 从 离散 相 到 连续 相 的 转变 , 使 
得 积累 的 拉 格 朗 日 液 膜 颗 粒 通过 VOF 的 形式 溶解 ， 质量、 动量 在 两 模型 之 间 传递 。 此 外 ， 离 散 粒子 撞击 辟 
膜 时 也 可 能 出 现 颗粒 反弹 和 液 膜 飞溅 现象 。 


(a) Before transition 


(b) During transition 


(c) After transition 


图 1 DTYV 转换 过 程 示意 图 
Fig.l Schematic Diagram of DTV Transition Process 
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2 ”模型 验证 和 边界 条 件 
2.1 DTYV 模型 验证 


在 液 滴 撞 击 壁面 的 数值 模拟 中 ，VOF 方法 能 够 捕捉 液 膜 的 细节 过 程 ， 但 对 网 格 精度 要 求 较 高 。 为 了 获 
得 一 定 程 度 的 计算 精度 ， 需 要 大 量 的 高 质量 网 格 ， 这 会 极 大 的 降低 计算 速度 。 相 比 之 下 ，DPM 模型 对 网 格 
精细 度 要 求 较 低 ， 拉 格 朗 日 粒子 可 以 被 放置 在 单个 大 网 格 中 ， 但 是 刚性 颗粒 不 便于 观察 碰撞 和 融合 细节 。 
因此 ， 本 研究 采用 了 基于 曲率 的 网 格 自 适应 细 化 方案 的 DTYV 方法 ， 以 平衡 计算 速度 和 计算 精度 ， 并 精确 
捕捉 相 界 面 的 演变 和 液 膜 铺展 细节 。 

为 验证 DTV 方法 的 可 靠 性 ， 采 用 了 两 液 滴 滴 落 融合 过 程 进行 验证 。 液 滴 材 质 为 87% 的 高 粘性 甘油 水 
液 滴 止 ， 直 径 3 mm， 密 度 1235 kg-m3。 这 些 液 滴 以 1.1 mys 的 速度 从 距离 壁面 3.4 mm 的 高 度 滴 落 。 液 滴 滴 
落 并 撞击 壁面 过 程 如 图 2 所 示 ， 在 撞击 壁面 之 前 ， 液 滴 一 直 以 DPM 离散 相 的 形式 存在 (a)， 撞 击 壁 面 形成 
拉 格 朗 日 壁 膜 (b)， 刚 性 DPM 颗粒 即使 撞击 壁面 也 不 会 发 生变 形 。 当 壁 膜 厚 度 超过 设 定 的 转换 值 时 ， 实 现 
了 由 DPM 到 VOF 相 的 转换 (c)， 同 时 质量 和 动量 发 生 传递 ， 使 得 VOF 液 滴 发 生 形 变 0(d)。 

图 2 还 展示 了 在 DPM 到 VOF 过 渡 过 程 中 的 网 格 细节 ， 基 于 曲率 的 网 格 自 适 应 细 化 多 次 触发 ， 形 成 了 
贴 合 VOF 液 面 的 新 的 气 液 边界 。 同 时 ， 在 气 液 界面 远离 之 后 ， 网 格 粗 化 也 会 触发 。 

将 两 液 滴 滴 落 融合 的 实验 图 像 和 模拟 过 程 进 行 比较 。 在 第 一 个 液 滴 滴 落 在 壁面 上 形成 稳定 形态 后 ， 将 
第 二 个 甘油 水 液 滴 滴 下 , 融合 和 发 展 过 程 如 图 3、4 所 示 , 它们 中 心 位 置 偏 移 量 分 别 为 0.73 mm 和 3.68 mm。 

当 离 散 相 颗粒 接触 VOF 液 相 后 并 达到 转换 条 件 后 ， 它 们 将 变 为 VOF 液 滴 并 继续 下 降 ， 与 第 一 个 液 滴 
发 生 相互 作 用 ， 直 到 合并 达到 稳定 状态 。 模 拟 结果 与 实验 图 片 BlG] 吻 合 较 好 ， 证 明了 DTV 模型 方法 的 可 行 
性 和 可 信 度 。 


(a) Before Impact (b) Before transition (c) After transition (d)Steady state 


2 液 滴 滴 落 撞击 壁面 的 转换 、 演 变 过 程 和 网 格 变化 情况 
Fig.2 Transitions, Evolutionary Processes, and Mesh Changes in Droplet Impacts on the Wall 
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Experiment Simulation Experiment Simulation 


0 ee. 


t=0 ms 
t=1 ms 


t=2 ms t=2 ms 


t=3 ms t=3 ms 


t=5 ms 


t=4 ms 


3 间距 0.73 mm 的 两 液 滴 滴 落 融合 过 程 BjD 对 比 4 间距 3.68 mm 的 两 液 滴 滴 落 融合 过 程 Blb 对 比 
Fig.3 Comparison of Droplet Falling and Fusion Fig.4 Comparison of Droplet Falling and Fusion 
Processes6B]6] with a Separation Distance of 0.73 mm ProcessesB}[5] with a Separation Distance of 3.68 mm 


2.2 ”计算 域 和 边界 条 件 


液 滴 射 流 撞 击 垂直 壁面 的 计算 域 如 图 4 所 示 , 其 中 高 度 肪 ”=50 mm 宽度 L ”=12 mm, 厚度 九 =3 mm 
利用 ICEM 对 计算 域 进行 了 结构 化 网 格 划 分 ， 并 对 壁面 进行 了 网 格 的 局 部 加 密 。 

计算 时 ， 选 择 VOF 模型 作为 两 相 流 模型 ，RNG k 一 gs 作为 汕 流 模型 。 选 择 single 注入 源 类 型 ， 注 射 
的 液 滴 为 液态 水 ， 直 径 为 0.7 mm， 密 度 为 998.2 kg:m”。 粒 子 速度 、 流 量 和 入 射 位 置 根据 不 同 工 况 进行 
设 定 ， 以 研究 不 同 射流 速度 、 不 同 射流 角度 对 壁面 液 膜 状 态 的 影响 规律 。 


第 一 作者 : 陈 青山 ， 男 ，1998 年 出 生 ， 博 士 研究 生 ， 核 反应 堆 热 工 水 力 
通讯 作者 : 田 文 喜 ，E-mail: wxtian@mail.xjtu.edu.cn 
收 稿 日 期 ，20XX-00-00， 修 回 日 期 ， 20XX-00-00 
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Droplet Injection 
不 


Wall 
1 


se 
D 
图 5 计算 域 几何 模型 


Fig.$ Geometric Representation of the Computational Domain 


3 ”计算 结果 分 析 


液 滴 以 不 同 速度 (流量) 和 角度 撞击 壁面 时 ， 形 成 的 液 膜 及 其 扩展 演变 情况 各 不 相同 。 本 文 分 析 了 注 
入 速度 和 入 射 角度 对 液 膜 扩 展演 变 过 程 的 影响 。 


3.1 射流 撞 壁 液 膜 特征 


连续 液 滴 以 DPM 颗粒 的 形式 、 以 一 定 的 角度 a 和 质量 流量 G 。 《人 入射 速度 wo ) 入 射 撞击 


垂直 壁面 ， 如 图 6， 撞击 点 为 0。DPM 粒子 在 接触 壁面 后 会 发 生 转 换 ， 形 成 VOF 液 膜 ， 而 后 续 颗 粒 在 接 
触 液体 后 也 会 发 生 转 换 。 由 于 撞击 导致 速度 方向 改变 ， 在 惯性 的 作用 下 ， 液 膜 从 撞击 点 O 向 四 周 扩 展 ， 其 


中 沿 着 壁面 向 下 的 质量 流量 为 C K 形成 汇流 。 此 外 ， 部 分 液 滴 在 撞击 壁面 后 会 发 生 溅 射 ， 脱离 壁面 形成 


小 液 滴 。 

液 滴 撞 击 壁 面 形成 液 膜 的 典型 铺展 情况 四 如 图 7 所 示 。 可 以 看 出 ， 液 滴 撞 击 壁 面 后 扩展 开 来 ; 由 于 能 量 
耗 散 和 壁面 粘性 作用 , 液 膜 动能 逐渐 减 小 ， 重 力作 用 增强 , 液 膜 向 下 扩展 ， 形 成 长 椭圆 形 的 液 膜 区 域 (PART 
D; 在 表面 张力 作用 下 , 液 膜 逐 渐 收 缩聚 拢 在 节点 K 处 ; 液 膜 汇聚 成 股 下 流 , 最 终 形 成 了 汇流 区 域 (PART ID。 
稳定 典型 液 膜 形状 类 似 “ 遇 虹 ”" 形 ， 其 中 HH 为 液 膜 区 高 度 ，W 为 液 膜 区 宽度 。 

液 膜 区 尺寸 是 影响 换 热 的 一 个 重要 因素 ， 研 究 中 通过 液 膜 区 高 度 和 宽度 来 表征 其 与 壁面 形成 的 接触 面 
面积 大 小 。 溅 射 率 也 是 高 速 液 滴 撞 击 壁面 时 的 又 一 重要 指标 ， 不 同 注入 速度 和 入 射 角度 都 会 引起 不 同 程度 
的 溅 射 ， 研 究 中 通过 比较 汇流 的 质量 流量 和 注入 液 滴 的 质量 流量 来 计算 获得 。 溅 射 率 y 的 计算 公式 外 为 ; 


G ye 
了 三 一 一 人 和 _x 1009%#(5) 
C 天 
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Before 
transition 


Sputtering : 
droplets 


After 
transition 


6 射流 撞 壁 示意 图 


7 典型 液 滴 射 流 撞 壁 实 验 效果 图 四 


Fig.6 Diagram of Jet Impinging on Wall Fig.7 Typical Experimental Image of Droplet Jet 


Impinging on WallI9] 


3.2 ”注入 速度 对 液 膜 特征 的 影响 分 析 
8 为 不 同 注入 速度 (45° 注 入 ， 注 入 速度 分 别 为 3 ms、S$m/s、8 m/s、10 m/s 和 12 m/s) 条 件 下 液 滴 


撞击 垂直 壁面 形成 的 液 膜 扩 展演 变 过 程 。 所 有 液 滴 具有 相同 的 撞击 距离 和 直径 ， 但 由 于 速度 不 同 ， 每 种 情 


况 下 液 滴 的 质量 流量 也 有 所 差异 。 液 滴 质量 流量 的 增加 直接 影响 了 液 膜 与 壁面 的 接触 面积 。 


注入 速度 对 液 膜 形 态 具 有 显著 影响 。 在 较 低 注入 速度 下 ， 形 成 的 液 膜 光滑 ， 汇 流 单一 有 旦 稳定 ; 随 着 注 


入 速度 的 增加 ， 液 膜 扩展 速度 加 快 ， 


液 接触 面积 增 大 。 速 度 增加 也 会 导致 液 膜 表面 扰动 增加 ， 形 成 更 多 的 


汇流 ， 使 液 膜 形态 发 生 较 大 变化 。 这 是 因为 高 速 液 滴 撞 击 壁 面 形成 的 液 膜 具有 更 高 的 初始 动能 ， 因 此 液 膜 
区 的 宽度 相 比 典型 情况 更 大 。 宽 度 的 增加 减弱 了 表面 张力 对 汇流 重力 的 影响 ， 导 致 液 膜 收缩 聚拢 的 节点 增 


多 ， 形 成 多 条 汇流 。 


图 9 显示 了 20 ms 时 的 液 膜 形态 和 与 壁面 的 接触 面 ， 随 着 速度 增加 ， 接 触 面积 、 液 膜 区 高 度 和 宽度 均 


明显 增加 ， 其 变化 趋势 如 图 10 所 示 


另 一 方面 ， 注 入 速度 的 增加 也 导致 液 滴 的 溅 射 率 增加 ， 变 化 曲线 如 图 10 所 示 。 较 大 的 注入 速度 不 仅 使 


液 滴 上 共有 更 大 的 向 四 周 扩 散 的 动能 ， 
滴 流 动 曲线 图 中 可 以 看 出 。 


o 


也 增加 了 垂直 于 壁面 方向 的 动能 ， 导 致 更 多 的 液 滴 溅 射 ， 从 图 11 的 液 


3.3 ”入 射 角度 对 液 膜 特征 的 影响 分 析 


图 12 为 不 同 入 射 角 度 (5 m/s 注 


E 入 ， 入 射 角度 分 别 为 30*>、45°、60°? 和 90°) 条 件 下 液 滴 撞 击 壁面 形成 


的 液 膜 扩展 演变 过 程 。 入 射 角度 越 小 ， 液 滴 在 撞击 壁面 后 紧 直 方向 的 速度 分 量 增 大 ， 导 致 汇流 流动 速度 增 
加 。 而 随 着 入 射 角度 的 增 大 ， 液 膜 速 度 的 竖 直 分 量 减 小 ， 同 时 水 平方 向 和 垂直 壁面 方向 速度 增 大 。 
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0ms lms 2ms 3ms 4ms Sms 6ms 7ms 8ms 9ms 10ms1Sms20ms 0ms lms 2ms 3ms 4ms Sms 6ms 7ms 8ms 9ms 10ms 15ms20ms 
(a) 3 m/s (b) 5 m/s 
、 4 
: . 
| 4 
0ms lms 2ms 3ms 4ms Sms 6ms 7ms 8ms 9ms lO0ms15ms20ms 0ms lms 2ms 3ms 4ms Sms 6ms 7ms 8ms 9ms 10ms15ms20ms 


c) 8 m/s i (d) 10 mys 


Oms lms 2ms 3ms 4ms Sms 6ms 7ms 8ms 9ms 10ms 15Sms20ms 


(e) 12 my/s 


图 8 不 同 注 入 速度 条 件 下 的 液 膜 特征 
Fig.8 Liquid Film Characteristics at Different Injection Rates 


(a)3my/s (b)Smy/s (c)8m/s (d)l0m/s (e)l2m/s 


图 9 不 同 注 入 速度 条 件 下 的 液 膜 形态 和 接触 面 示意 图 
Fig.9 Schematic Diagram of Liquid Film Morphology and Contact Surfaces at Different Injection Rates 
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3.1m/s 5.1mAs 8.lm/s 10.1nvs 12.1mys 
六 
E 避 
E E 
己 er 
三 
品 号 
a 
只 
Qnys ‘Qnys Onvs Onms Onvs 
Droplets velocity /ms (a)3m/s (b)Sm/s (c)8m/s (d)10m/s (e)l2m/s 
图 10 不 同 注 入 速度 条 件 下 的 液 膜 特征 曲线 图 11 不 同 注 入 速度 条 件 下 的 液 滴 流 动 轨迹 示意 图 
Fig.10 Characteristic Curve of Liquid Film at Different Fig.11 Schematic Diagram of Droplet Flow Trajectory at 
Injection Rates Different Injection Rates 


由 图 13 可 见 ， 液 膜 速度 分 量 的 变化 直接 影响 了 液 膜 表 面 的 稳定 性 ， 使 得 液 膜 形态 不 再 规则 ， 并 且 液 膜 
与 壁面 的 接触 面积 也 发 生变 化 。 液 膜 区 宽度 随 着 入 射 角度 的 增加 而 增 大 ， 而 高 度 则 随 着 入 射 角度 的 增 大 而 
减 小 ， 如 图 14 所 示 ， 宽 度 的 增加 速率 略 高 于 高 度 的 减 小 速率 。 

液 滴 溅 射 率 也 随 着 垂直 于 壁面 的 速度 分 量 的 增加 而 增加 ， 如 图 15 所 示 。 在 较 小 的 入 射 角度 区 间 ， 溅 射 
率 的 增加 相对 较 小 ， 而 在 较 大 的 区 间 ， 增 加 速率 更 快 。 特 别 是 在 90” 入 射 的 情况 下 ， 液 滴 垂 直 撞 击 紧 直 壁 
面 ， 导 致 73.65% 的 液 滴 发 生 反 弹 、 接 触 壁 面 后 又 脱离 ， 使 得 液 滴 溅 射 率 远 高 于 其 他 工 况 。 


3 7 20 


Oms lms 2ms 3ms 4ms Sms 6ms 7ms 8ms 9ms 10ms 1Sms 20ms Oms lms 2ms 3ms 4ms Sms 6ms 7ms 8ms 9ms 10ms15ms20ms 

(a) 

, 4 
, | 
le 
Oms lms 2ms 3ms 4ms Sms 6ms 7ms 8ms 9ms 10ms15ms20ms 0ms lms 2ms 3ms 4ms Sms 6ms 7ms 8ms 9ms 10ms 1Sms 20ms 
O o 
(c) 60 (d) 90 


图 12 不 同 入 射 角度 条 件 下 的 液 膜 特征 
Fig.12 Liquid Film Characteristics at Different Incidence Angles 
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(a)30° (b)45° (c)60° (d)90° 


图 13 不 同 入 射 角度 条 件 下 的 液 膜 形态 和 接触 面 示意 图 
Fig.13 Schematic Diagram of Liquid Film Morphology and Contact Surfaces at Different Incidence Angles 


2 ee 
10 上 
名 上 
入 
| Jo 
2 
30 40 50 60 70 80 90 ji 


1 
ES 


Length / mm 
大 
Sputtcring rate / % 


La 
So 


Injection angle /° 全 人 
图 14 不 同 入 射 角度 条 件 下 的 液 膜 特征 曲线 图 15 不 同 入 射 角度 条 件 下 的 液 滴 流 动 轨 迹 示意 图 
Fig.14 Characteristic Curve of Liquid Film at Different Fig.1S Schematic Diagram of Droplet Flow Trajectory at 
Incidence Angles Different Incidence Angles 


不 同 注 入 速度 和 不 同 入 射 角度 对 液 膜 特征 的 影响 规律 与 前 人 研究 结果 吻合 较 好 。 此 外 ， 相 比 于 其 他 的 
耦合 方法 ， 本 文 研 究 中 采用 的 DTV 方法 在 确保 计算 精度 的 前 提 下 ， 表 现 出 更 快 的 计算 速度 ， 能 够 有 效 地 
实现 离散 相 到 连续 相 的 耦合 。 因 此 ， 该 方法 可 用 于 分 析 和 研究 液 滴 射流 对 壁面 、 液 面 以 及 壁 膜 的 撞击 情况 


4 ”结论 


本 文采 用 DTYV 方法 对 液 滴 撞击 垂直 壁面 后 的 液 膜 扩展 和 演变 过 程 进 行 了 分 析 ， 并 研究 了 注入 速度 和 
入 射 角度 对 液 膜 形态 和 液 滴 溅 射 率 的 影响 规律 。 得 到 的 主要 结论 如 下 : 

(1)DTYV 模型 适用 于 离散 液 滴 撞击 壁面 、 液 面 、 液 膜 过 程 的 模拟 ,能 够 准确 计算 撞击 过 程 中 液 膜 融合 和 
液 膜 演变 。 

(2) 增 大 注入 速度 有 利于 液 膜 的 扩展 、 增 大 液 膜 与 壁面 的 接触 面积 、 提 高 液 膜 区 高 度 和 宽度 ， 但 也 会 增 
加 液 滴 溅 射 率 。 

(3) 入 射 角度 对 液 膜 形态 有 显著 影响 ， 增 大 入 射 角度 能 够 增 大 液 膜 区 宽度 、 减 小 液 膜 高 度 ， 同 时 导致 液 
滴 溅 射 率 变化 速度 加 快 ， 特 别 是 在 90" 时 ，73.65% 的 液 滴 发 生 溅 射 脱离 壁面 。 
在 安全 壳 、 稳 压 器 和 其 他 核 动力 设备 表面 的 喷 淋 降温 过 程 中 ， 增 大 液 膜 与 壁面 的 接触 面积 是 强化 降温 
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F DPM-to-VOF 的 液 滴 射 流 撞 壁 过 程 研究 


的 一 个 手段 。 液 体 注 入 速度 和 入 射 


传 热 模型 进行 喷 淋 状态 下 壁面 降温 下 
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